When the French mathematician Joseph Fourier (1768-1830) was trying to solve a prob¬ 
lem in heat conduction, he needed to express a function / as an infinite series of sine and 
cosine functions: 


Fourier Series 


m 


fix) — ao + 2 i a n cos nx + b n sin nx) 


n=l 


= ao + a\ cos x + ci 2 cos 2x + a 3 cos 3x + 


+ b\ sin x + b 2 sin 2x + b 3 sin 3x + ■ ■ ■ 


Earlier, Daniel Bernoulli and Leonard Euler had used such series while investigating prob¬ 
lems concerning vibrating strings and astronomy. 

The series in Equation 1 is called a trigonometric series or Fourier series and it turns 
out that expressing a function as a Fourier series is sometimes more advantageous than 
expanding it as a power series. In particular, astronomical phenomena are usually periodic, 
as are heartbeats, tides, and vibrating strings, so it makes sense to express them in terms 
of periodic functions. 

We start by assuming that the trigonometric series converges and has a continuous func¬ 
tion fix) as its sum on the interval [—77, 77], that is. 


00 


fix) = ao + 2 i a n cos nx + b„ sin nx) 


— 7T =£ X *£ 77 


Our aim is to find formulas for the coefficients a„ and b„ in terms of /. Recall that for a 
power series fix) = X c„(x — a)' 1 we found a formula for the coefficients in terms of deriv¬ 
atives: c„ = f (n \a)/n\. Here we use integrals. 

If we integrate both sides of Equation 2 and assume that it’s permissible to integrate the 
series term-by-term, we get 




But 



1 

cos nx dx = — sin nx 
n 


n 


1 


[sin/777 — sin(— H77)] = 0 


because n is an integer. Similarly, sin nx dx = 0. So 



1 
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llll Notice that a a is the average value of / over 
the interval \-tt, it}. 


and solving for a 0 gives 



To determine a„ for « 5* 1 we multiply both sides of Equation 2 by cos mx (where m is 
an integer and m S 5 1) and integrate term-by-term from —77 to 77 : 


a 0 + 2 (fln cos nx + b„ sin nx) 

n= 1 


cos mx dx 


j f{x) cos mx dx = j" 

|~4~| = flo cos mx dx + 2 cos nx cos mx dx + 2 sin nx cos mx dx 

J-ir n=l •' _,r /i=l 

We’ve seen that the first integral is 0. With the help of Formulas 81, 80, and 64 in the Table 
of Integrals, it’s not hard to show that 

I sin nx cos mx dx = 0 for all n and m 

J —7T 

0 for n 7 ^ m 


f 

J —IT 


cos nx cos mx dx = 


77 for n = m 

So the only nonzero term in (4) is a„,Tr and we get 

j /(x) cos mx dx = a m 77 
Solving for a m , and then replacing m by n, we have 


H 


we get 




a n = — 

77 J 

/(x) cos nx dx n = 1, 2, 3, ... 

77 

we multiply both sides of Equation 2 by sin mx and integrate from 

b n = - 

77 

/(x) sin nx dx n = 1, 2, 3,... 

— 77 


We have derived Formulas 3, 5, and 6 assuming / is a continuous function such that 
Equation 2 holds and for which the term-by-term integration is legitimate. But we can still 
consider the Fourier series of a wider class of functions: A piecewise continuous function 
on [a, b] is continuous except perhaps for a finite number of removable or jump disconti¬ 
nuities. (In other words, the function has no infinite discontinuities. See Section 2.5 for a 
discussion of the different types of discontinuities.) 
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llll Engineers use the square-wave function in 
describing forces acting on a mechanical system 
and electromotive forces in an electric circuit 
(when a switch is turned on and off repeatedly). 
Strictly speaking, the graph of / is as shown 
in Figure 1(a), but it's often represented as in 
Figure 1(b), where you can see why it's called a 
square wave. 


[71 Definition Let / be a piecewise continuous function on [—77, 77]. Then the 
Fourier series of / is the series 


a 0 


+ 2 ( a n cos nx + b„ sin nx) 


where the coefficients a„ and b„ in this series are defined by 


ao=— r fix) 

ATT J-n 


dx 


1 j* 57 

a„ = — fix ) cos nx dx 

77 J-TT 


b„ = — | * fix) si 

77 J-TT 


sin nx dx 


and are called the Fourier coefficients of /. 


Notice in Definition 7 that we are not saying fix) is equal to its Fourier series. Later we 
will discuss conditions under which that is actually true. For now we are just saying that 
associated with any piecewise continuous function / on [—77, 77] is a certain series called 
a Fourier series. 

EXAMPLE 1 Find the Fourier coefficients and Fourier series of the square-wave function 
/ defined by 


, . Jo if — 77 ^ X < 0 , . . . 

fix) = j and fix + 2 77) = fix) 

1 if 0 x < 77 


So / is periodic with period 2tt and its graph is shown in Figure 1. 
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(b) 


SOLUTION Using the formulas for the Fourier coefficients in Definition 7, we have 


ao — 


1 fir , , 1 70 1 fir 1 , . 

fix) dx = — 0 dx + — 1 dx = 0 + — (77) 

277 J-v 277 J-n 277 •'<> 277 
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FIGURE 1 

Square-wave function 
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and, for n 3= 1, 


1 fir 1 ro 1 ff 

— /(x) cos nx dx = — 0 dx H- 

77 «/ — 7 r 77 J - 77 77 JO 


cos nx dx 


= 0 + 


1 


sin nx 



77 n 


-(sin 72-77 — sin 0) = 0 

mr 


llll Note that cos ntr equals 1 if n is even 
and -1 if n is odd. 


bn = — r f{x) si 

TT J-TT 


sin nx dx = — 


1 ro 1 r-n 

— 0 dx H-si 

77 J-TT 77 JO 


sin x dx 


1 cos nx 


-(COS 77-77 — COS 0) 

n 77 


0 if n is even 
2 

- if n is odd 

mr 


The Fourier series of / is therefore 


ao + a\ cos x + «2 cos 2x + a 3 cos 3x + • • ■ 


+ b i sin x + bi sin 2x + i> 3 sin 3x + 


1 

= — + 0 + 0 + 0 + -" 
2 


H-sin x + 0 sin 2x H-sin 3x + 0 sin 4x H-sin 5x + 

77 377 577 


1 2 2 2 2 

=-1-sin x H-sin 3x H-sin 5x H-sin lx + 

2 77 377 577 777 


Since odd integers can be written as n = 2k — 1, where k is an integer, we can write the 
Fourier series in sigma notation as 


2 


+ 2 


2 

(2k — 1)77 


sin(2£ — l)x 


In Example 1 we found the Fourier series of the square-wave function, but we don’t 
know yet whether this function is equal to its Fourier series. Let’s investigate this question 
graphically. Figure 2 shows the graphs of some of the partial sums 


S„(x) 


12 2 2 

=-1-sin x H-sin 3x + ■ • • H-sin nx 

2 77 377 n 7T 


when n is odd, together with the graph of the square-wave function. 
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FIGURE 2 Partial sums of the Fourier series for the square-wave function 

We see that, as n increases, S„(x) becomes a better approximation to the square-wave 
function. It appears that the graph of S„(x) is approaching the graph of fix), except where 
x = 0 or x is an integer multiple of 77. In other words, it looks as if / is equal to the sum 
of its Fourier series except at the points where / is discontinuous. 

The following theorem, which we state without proof, says that this is typical of the 
Fourier series of piecewise continuous functions. Recall that a piecewise continuous func¬ 
tion has only a finite number of jump discontinuities on [—77, 77]. At a number a where / 
has a jump discontinuity, the one-sided limits exist and we use the notation 

f(a + ) = lim fix) f(a~) = lim fix) 


| 1 T] Fourier Convergence Theorem If / is a periodic function with period 2i 7 and / and 
f are piecewise continuous on [ — tt, tt ], then the Fourier series (7) is convergent. 
The sum of the Fourier series is equal to f(x) at all numbers x where / is continu¬ 
ous. At the numbers x where / is discontinuous, the sum of the Fourier series is 
the average of the right and left limits, that is 

\[f[x + ) +/(*“)] 


If we apply the Fourier Convergence Theorem to the square-wave function / in 
Example 1, we get what we guessed from the graphs. Observe that 

/(0 + ) = lim fix) = 1 and /(0~) = lim fix) = 0 

x— »0 + x —>0 

and similarly for the other points at which / is discontinuous. The average of these left and 
right limits is so for any integer n the Fourier Convergence Theorem says that 


2 


+ i 


2 

(2 k — 1)77 


sin(2k — 1 )x = 



if n ¥= mr 

if X = HTT 


(Of course, this equation is obvious for x = iitt.) 
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Functions will Period PL 


If a function / has period other than 2i 7 , we can find its Fourier series by making a change 
of variable. Suppose fix) has period 2L, that is fix + 2 L) = fix) for all x. If we let 
t = 77 x/L and 


git) = fix) = fi.Lt/ it) 

then, as you can verify, g has period 2tt and x = ±L corresponds to t — ±ir. The Fourier 
series of g is 

a 0 + 2 i a n cos nt + b„ sin nt) 

n= 1 

where 

a 0 = r 3(f) dt 

Lit d-TT 


a 


n 



cos nt dt 


b„ 



sin nt dt 


If we now use the Substitution Rule with x = L f it, then t = ttx/L, dt = ( tt/L) dx, and 
we have the following 


llll Notice that when L = v these formulas are 
the same as those in (7). 


[~9~1 If / is a piecewise continuous function on [— L, L\, its Fourier series is 


ao + 2 

n=i 


a „ cos 


+ b„ sin 


where 


and, for n S 5 1, 


ao = - f fix) dx 

2L J-L J v 


1 fL 


&n l 


L m 


COS 


n 7TX 


dx 


1 ^ 
b " = ~L J- 


j fix) sin 


n 77 v 


dx 


Of course, the Fourier Convergence Theorem (8) is also valid for functions with period 
2 L. 

EXAMPLE 2 Find the Fourier series of the triangular wave function defined by fix) = \ x \ 
for —1 ^ x 1 and fix + 2) = fix) for all x. (The graph of / is shown in Figure 3.) 
For which values of x is fix) equal to the sum of its Fourier series? 



FIGURE 3 

The triangular wave function 
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llll Notice that ao is more easily calculated as 
an area. 


SOLUTION We find the Fourier coefficients by putting L = 1 in (9): 


ao = a j" \x\dx = ^j (~x)dx + \ j" x dx 


= — -r 2 l° + = - 

4-* J_[ T 4 X jo 2 


and for n 5= 1, 


Qn 


j" | x | cos(n7rx) dx = 2 j x cos(h77x) dx 


because y = \ x \ cos(/nrx) is an even function. Here we integrate by parts with u = x 
and dv = cos(n 7rx) dx. Thus, 


a„ = 2 


-sin(«77x) 

H7T 


= 0 - 


n 77 


COs(u 77x) 


1VTT 

1 


f sin(«7rx) 
Jo 


dx 


nir 


-(cos tl 77 — 1) 


0 


Since y = |x| sin(«7rx) is an odd function, we see that 

b n = J \x \ sin(«7rx) dx = 0 

We could therefore write the series as 

1 ^ 2(cos HIT — 1) 

-h 2j - 2 — 2 -COS(« TTX) 

2 „=1 n it 

But cos mr = 1 if n is even and cos nir = — 1 if n is odd, so 


a n = — —x(cos iut — 1) = “ 

IV 7T 


0 if n is even 

4 

-r—- if n is odd 


Therefore, the Fourier series is 


1 4 , N 4 4 , 

-xcos(77\r)- t cos(3 77x)- — COs(5 77X') 

2 7 t 9v- 25 tt~ 




2 „=1 (2k - 1)V 


-cos((2£ — 1) 7 tx) 


The triangular wave function is continuous everywhere and so, according to the Fourier 
Convergence Theorem, we have 

fix) = - S —— 4 , , cos((2 k - 1)t tx) for all x 

2 „=i (2k - 1) 7r 
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In particular, 

1 » 4 

\x\ = — - S ~r~ - —r~rcos((2k - 1)t tx) for -1 =S x =S 1 

2 k = 1 (2k - 1 )~7T~ 


Fourier Series and Music 


One of the main uses of Fourier series is in solving some of the differential equations that 
arise in mathematical physics, such as the wave equation and the heat equation. (This is 
covered in more advanced courses.) Here we explain briefly how Fourier series play a role 
in the analysis and synthesis of musical sounds. 

We hear a sound when our eardrums vibrate because of variations in air pressure. If a 
guitar string is plucked, or a bow is drawn across a violin string, or a piano string is struck, 
the string starts to vibrate. These vibrations are amplified and transmitted to the air. The 
resulting air pressure fluctuations arrive at our eardrums and are converted into electrical 
impulses that are processed by the brain. How is it, then, that we can distinguish between 
a note of a given pitch produced by two different musical instruments? 

The graphs in Figure 4 show these fluctuations (deviations from average air pressure) 
for a flute and a violin playing the same sustained note D (294 vibrations per second) as 
functions of time. Such graphs are called waveforms and we see that the variations in air 
pressure are quite different from each other. In particular, the violin waveform is more 
complex than that of the flute. 


FIGURE 4 

Waveforms 




We gain insight into the differences between waveforms if we express them as sums of 
Fourier series: 


7 it 


P(t) = ao + ai cos I — I + bi sin I — I + ai cos 


TTt 


2tt t 


+ £>2 sin 


2irt 


+ 


In doing so, we are expressing the sound as a sum of simple pure sounds. The difference 
in sounds between two instruments can be attributed to the relative sizes of the Fourier 
coefficients of the respective waveforms. 

The nth term of the Fourier series, that is. 


a„ cos 


mrt 


+ b,, 


mrt 


is called the /zth harmonic of P. The amplitude of the nth harmonic is 


A n = \Jal + bl 

and its square, A 2 „ = al + bl, is sometimes called energy of the nth harmonic. (Notice that 
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for a Fourier series with only sine terms, as in Example 1, the amplitude is A„ = | b„ | and 
the energy is A\ = b;,.) The graph of the sequence {,4;,} is called the energy spectrum of 
P and shows at a glance the relative sizes of the harmonics. 

Figure 5 shows the energy spectra for the flute and violin waveforms in Figure 4. Notice 
that, for the flute, Ar„ tends to diminish rapidly as n increases whereas, for the violin, the 
higher harmonics are fairly strong. This accounts for the relative simplicity of the flute 
waveform in Figure 4 and the fact that the flute produces relatively pure sounds when 
compared with the more complex violin tones. 
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FIGURE 5 

Energy spectra 

(a) Flute 

(b) Violin 


In addition to analyzing the sounds of conventional musical instruments, Fourier series 
enable us to synthesize sounds. The idea behind music synthesizers is that we can combine 
various pure tones (harmonics) to create a richer sound through emphasizing certain 
harmonics by assigning larger Fourier coefficients (and therefore higher corresponding 
energies). 


Exercises 


1-6 nil A function / is given on the interval [—77, 7r] and / is 
periodic with period 277. 

(a) Find the Fourier coefficients of /. 

(b) Find the Fourier series of /. For what values of x is fix) equal 
to its Fourier series? 

@3 ( c ) Graph / and the partial sums S 2 , S 4 , and S 6 of the Fourier series. 


1. 

2 . 

3. 

4. 

5. 


fix) 

fix) 


r 1 t — 77 =s * < 0 

[ — 1 ifO®Sx<T7 

(" 0 if — 77 =£ x < 0 
lx ifO=Sx<77 


fix) = x 
fix) = X 1 


fix) 


To if — 77 =S x < 0 
[cos x if 0 =£ x < 77 


— 1 if — 77 =£ x < — 77/2 
1 if — 77/2 =S x < 0 

0 if 0 =£ x < 77 


7-11 llll Find the Fourier series of the function. 
1 if I x I < 1 


7. fix) = 


8 - fix) = 


9. f(x) = 

■' [0 if 0 « x < 4 

10 . f(x) = 1 - x, -1 =£ x =s 1 

11. fit) = sin(377f), — 1 =S t =£ 1 


0 if 1 =£ | x | < 2 

0 if -2 x < 0 
1 if 0 =£ x < 1 
0 if I =S x < 2 

-x if —4 x < 0 


fix + 4) = fix) 

fix + 4) = fix) 

fix + 8) = fix) 
fix + 2) = fix) 


12. A voltage E sin cot, where t represents time, is passed through a 
so-called half-wave rectifier that clips the negative part of the 
wave. Find the Fourier series of the resulting periodic function 


fit) = 


0 


if S£ f < 0 

O) 


E sin cot if 0 ^ < — 


6- fix) 


fit + Itt/co) = fit) 
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13-16 llll Sketch the graph of the sum of the Fourier series of / 
without actually calculating the Fourier series. 


13. f(x) 


f-1 if —4 =£ jc < 0 
[3 if 0 =S x < 4 


14. f{x) 


fx if — 1 X < 0 

]_ 1 — x if 0 jc < 1 


15. f{x) = x\ -l=Sx=sl 

16. fix) = e\ -2«x=s2 


17. (a) Show that, if — 1 =£ jc 1, then 


x 


2 


— + 2 ( — 1)" ~~i rcos(n7TJr) 
3 „_i n-n~ 


(b) By substituting a specific value of x, show that 


i 


1 


n 


2 


77 


2 


6 


18. Use the result of Example 2 to show that 


1 



IT 


2 


8 


19. Use the result of Example 1 to show that 


1 - 



1 

-j_ . . . 

7 


7 T 

4 


20. Use the given graph of / and Simpson’s Rule with n = 8 to 
estimate the Fourier coefficients a 0 , a u cii, b 1 , and b 2 . Then use 
them to graph the second partial sum of the Fourier series and 
compare with the graph of /. 





















